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NOMENCLATURE 


Preexponential  constant  in  Arrhenius  rate  of  solid-phase 
reaction. 

Preexponential  constant  in  Arrhenius  rate  of  gas-phase 
reaction. 

Solid  heat  capacity  per  unit  mass. 

Gas  specific  heat  at  constant  pressure. 

Diffusion  coefficient. 

Activation  energy  for  gas-phase  reaction. 

Activation  energy  for  solid-phase  reaction. 

Heat  released  in  solid  interior  by  main  decomposition  reaction, 
per  unit  mass  of  reactant  consumed. 

Heat  released  in  solid  interior  by  stabilizer  reaction,  per  unit 
mass  of  reactant  consumed. 

Heat  released  in  gas-phase  reaction  per  unit  mass  of  reactant 
consumed. 

Parameter  defined  by  Eq  (28) . 

Sum  of  heat  of  gasification  and  heat  conducted  into  interior  of 
solid,  both  per  unit  mass  of  solid  gasified. 

Heat  of  gasification  per  unit  mass  of  solid  at  temperature  Tg. 

Mass  burning  rate. 

Pressure  exponent  of  burning  rate. 

Steady  combustion  pressure. 

External  effective  radiant  energy  flux. 

External  radiant  energy  flux  absorbed  in-depth  by  solid. 

External  radiant  energy  flux  absorbed  by  solid  gasification  layer. 
Universal  gas  constant. 

Regression  rate  of  solid. 

Temperature. 

Flame  temperature. 

Initial  temperature  of  solid. 

Temperature  of  solid  surface. 

Time. 

Distance  normal  to  solid  surface. 

Distance  between  solid  surface  and  gas-phase  reaction  zone. 


NOMENCLATURE  (Cont'd) 


Y Reactant  mass  fraction, 

a Gas  diffusivity, 

A Burning  rate  eigenvalue, 

A Solid  thermal  conductivity. 

A Gas  thermal  conductivity. 

o 

y Absorption  coefficient  for  solid, 

v Reaction  order  of  gas-phase  reaction. 

5 Nondimensional  distance  in  solid. 

JIT  Temperature  sensitivity  coefficient  of  burning  rate, 

p Solid  density, 

p Gas  density. 

o 

a Nondimensional  absorption  coefficient  for  solid. 

t Nondimensional  solid  temperature. 

Subscript 


o 


Denotes  absence  of  external  radiant  energy  flux 


I.  INTRODUCTION 


A new  theoretical  analysis  aimed  at  developing  a better  understanding 
of  the  influence  of  thermal  effects  in  solid  propellant  combustion  is 
given  here.  It  also  aims  at  obtaining  burning  rate  predictions  which 
exhibit  explicitly  the  dependence  on  radiant  flux.  Practical  motivation 
ranges  from  interest  in  burning  rate  control  to  concern  about 
detonation  suppression. 

Solid  propellant  burning  is  sustained  by  several  modes  of  energy 
release  and  transfer.  Meat  generation  in  the  solid  phase  and  heat 
feedback  from  the  gas-phase  flame  to  the  burning  surface  are  major 
driving  mechanisms,  depending  upon  which  process  is  controlling.  The 
study  of  thermal  effects,  derived  by  stimulating  the  burning  propellant 
with  an  external  source  of  radiant  energy,  constitutes  a fundamental 
method  of  testing  the  applicability  of  combustion  models,  and  lias  the 
experimental  advantage  of  readily  available  equipment.  This  method  also 
has  the  theoretical  advantage  that  many  important  interpretations  can  be 
made  without  having  to  resort  to  uncertain  descriptions  of  flame  structure 
or  schemes  about  the  intricate  chemical  reactions  taking  place.  The 
analysis  is  based  on  energy  considerations  with  no  direct  inclusion  of 
chemistry;  i.e.,  the  chemical  reactions  involved  are  assumed  to  proceed 
faster  with  added  thermal  energy  but  not  differently. 

Because  of  the  more  immediate  interest  of  double-base  propellants 
in  Army  applications,  this  type  of  solid  propellant  has  been  investigated 
first.  Composite  propellants  may  be  considered  later  in  the  study  and 
have  not  been  included  here.  Therefore,  the  multiplicity  of  conditions 
that  characterize  the  combustion  of  double-base  propellants  within 
rocket  operating  pressures  have  been  examined  in  detail.  The  descriptions 
of  the  interaction  between  the  thermal  radiation  and  the  combustion  zone 
result  from  limiting  procedures  involving  the  three  principal  solid- 
phase  lengths  of  surface  reaction,  heat  conduction,  and  radiation 
absorption.  By  allowing  the  absorption  length  to  vary  over  a wide  range, 
the  effect  of  external  radiant  energy  is  determined  for  each  condition 
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of  interest.  When  the  absorption  length  becomes  very  small,  surface 
gasification  appears  to  dominate  and  no  effect  from  the  external 
radiation  is  predicted,  provided  that  the  gas-phase  flame  is  not 
controlling.  Generation  of  heat  in  the  interior  of  the  solid  phase  can 
also  become  an  important  factor  when  the  radiation  is  absorbed  deep 
inside  the  solid. 

To  simplify  the  presentation  of  results,  the  thermal  stimulus 
provided  by  the  external  source  is  represented  by  the  net  extra  flux 
which  would  be  effectively  absorbed  at  the  burning  propellant  surface 
after  allowing  for  losses  along  the  optical  path.  This  indirectly 
supposes  that  the  product  gases  are  transparent,  which  is  of  course  not 
true.  Metallized  propellants  are  excluded  from  the  experiment  for  this 
reason,  and  not  because  of  any  theoretical  requirement.  The  surface 
gasification  and  gas  flame  reaction  are  represented  by  Arrhenius  rate 
expressions  assumed  to  have  a high  activation  energy.  The  surface  and 
final  flame  temperatures  are  viewed  in  the  traditional  sense  as  free 
parameters  arrived  at  by  energy  balance.  No  assumptions,  however,  are 
introduced  regarding  the  temperature  profile  between  the  surface  and 
final  temperatures  or  regarding  other  elements  of  flame  structure. 

In  the  past,  the  consideration  of  thermal  effects  in  solid  propellant 

combustion  evolved  from  a very  specialized  practical  problem.  It 

concerned  the  degree  of  augmentation  of  burning  rate  in  rocket  chambers 

1* 

as  compared  to  the  strand-burner  rate.  The  augmentation  observed  for 
the  same  pressure  was  explained  in  terms  of  the  thermal  radiation  field 
generated  inside  the  rocket  chamber.  From  this  early  work  comes  the 
hypothesis  of  equivalence  between  heat  absorbed  and  increase  in  initial 
propellant  temperature.  The  present  study  has  confirmed  the  appropriateness 
of  this  "equivalence  principle",  and  has  extended  its  validity  through 
a rigorous  treatment  of  the  problem.  It  has  further  yielded  the  very 
important  theoretical  result  of  delineating  the  unique  burning  conditions 
under  which  the  equivalence  principle  fails  to  apply. 

* References  are  listed  on  page  55. 
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Another  practical  aspect  of  the  thermal  effect  is  the  temperature 
sensitivity  of  solid  propellant  combustion.  As  a general  problem,  the 
temperature  dependence  is  a primary  propellant  parameter  (likewise  the 
pressure  dependence) , which  cannot  be  defined  without  highly  accurate 
combustion  models,  currently  nonexistant.  Conceptually,  however,  this 
problem  is  a part  of  the  general  theory  of  thermal  effects.  The  findings 
of  this  work  are  also  significant  in  this  regard,  by  showing  the  strong 
effect  that  solid  phase  heat  generation  from  individual  reactants,  like 
stabilizers,  can  have  on  temperature  sensitivity. 

The  main  result  can  be  summarized  as  follows.  Somewhere  in  the 
overall  combustion  zone  there  is  a kinetic  process  which  determines  the 
burning  rate;  the  radiation  absorbed  after  this  control  zone  is 
equivalent  to  an  increase  in  initial  propellant  temperature;  the 
radiation  absorbed  outside  is  ineffective.  The  exception  to  this  rule 
is  the  case  of  absorption  occurring  right  at  the  surface  reaction  layer 
accompanied  by  a change  in  surface  temperature,  and  provided  that  the 
gas-phase  flame  is  not  controlling  the  burning  rate.  Extremely  opaque 
solid  propellants  burning  at  low  pressures  probably  fall  in  this  category. 
Also,  if  the  external  thermal  radiation  were  to  cause  substantial 
photochemical  changes  in  the  combustion  process,  a different  type  of 
exception  would  appear.  No  positive  evidence  has  been  found  for  this 
phenomenon,  however.  As  for  the  experimental  verification  of  results, 
double-base  propellants  burning  in  a linear  regime  at  the  usual  rocket 
pressures  are  the  logical  choice.  Nonlinear  burning  propellants,  i.e., 
propellants  exhibiting  pronounced  mesa  or  plateau  regions,  have 
complicated  and  unknown  pressure  and  temperature  dependences  that  need 
to  be  resolved  first  by  appropriate  combustion  modelling. 

II.  GAS-PHASE  FLAME  REACTION  AND  SOLID-PHASE  RESPONSE 
A.  Deposition  of  External  Radiation 

The  analysis  is  based  on  the  one-dimensional  burning  of  a 
homogeneous  solid  propellant.  The  solid  phase  occupies  the  region  x<0 
and  the  gas-phase  the  region  x>0;  external  radiation  is  incident 
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on  the  propellant  from  the  gas.  Because  of  losses  along  the  optical 
path,  some  attenuation  of  the  initial  radiation  intensity  occurs. 

Losses  by  gas  absorption  can  be  significant,  especially  in  metallized 
propellants,  as  well  as  losses  by  reflection  at  the  propellant  surface. 
Emission  of  radiation  by  the  solid  will  usually  be  a small  correction. 

In  principle  the  analysis  is  applicable  in  the  presence  of  losses; 
it  is  only  necessary  to  define  the  radiant  flux,  Q,  as  the  one  that 
actually  reaches  the  propellant  surface. 

Transport  of  radiation  within  the  solid  phase  is  a truly  complicated 
phenomenon.  Deposition  of  radiation  by  absorption  and  scattering  are 
very  important  for  real  propellants;  however,  a detailed  optical  and 
mathematical  treatment  would  remain  far  from  practical  application. 
Therefore,  both  processes  are  lumped  into  an  empirical  Beer's  law 

liX 

relationship,  Qe  , with  an  effective  overall  absorption  coefficient  p, 
to  be  determined  experimentally.  Scattered  radiation  leaving  the  surface 
is  counted  as  reflection.  This  assumed  law  has  the  right  general 
functional  form  and  thus  provides  a good  practical  description  of  the 
phenomenom. 

Other  assumptions  pertain  to  the  physical  properties.  The 

2 

dependence  of  p on  the  wavelength  of  the  radiation  is  not  simple  , but 
it  is  preferable  to  use  a value  of  p corresponding  to  the  wavelength 
range  of  most  of  the  radiant  flux.  There  should  be  no  difficulty  in 
doing  this  unless,  for  example,  half  of  the  radiation  is  absorbed  deep 
in-depth  and  the  other  half  at  the  surface.  Anyway,  p is  taken  to  be 
a known  constant.  Experimentally  determined  values  of  p range  from 
about  10  cm"*  for  fairly  transparent  to  almost  200  cm"*  for  very  opaque 
propellants,  80  to  100  cm*  being  more  often  the  average  for  double-base 
propellants.  For  simplicity,  the  heat  capacity,  thermal  conductivity, 
and  density  of  the  solid  are  assumed  constant.  Corresponding  assumptions 
are  made  for  the  gas,  but  these  are  stated  in  Appendix  A,  giving  the 
analysis  of  the  gas-phase  reaction  zone.  The  solid  phase  is  considered 
first,  since  that  is  where  the  radiation  makes  itself  felt  initially. 
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B.  Solid-Phase  Temperature  Profiles 

The  treatment  of  the  thermal  response  of  the  solid  phase  without 
in-depth  heat  generation  includes  the  further  assumption  that  the  length 
of  the  surface  layer,  where  chemical  reaction  occurs,  is  very  small 
compared  with  either  the  lengths  of  the  heat  conduction  zone  or  the 
radiation  absorption  zone.  If  gasification  is  a true  interface  process 
(e.g.,  vaporization),  then  there  is  no  approximation  at  all.  However, 
although  propellants  have  a polymeric  structure,  and  polymers  tend  to 
gasify  by  in-depth  degradation,  it  has  been  shown  that  even  in  this  case, 
and  because  of  the  high  activation  energy  for  degradation,  the  reaction 
zone  length  is  typically  less  than  one-tenth  the  length  of  the  heat 
conduction  zone.  Except  for  propellants  of  unusually  high  opacity,  it 
is  unlikely  that  the  radiation  can  be  absorbed  strongly  enough  for  most 
of  it  to  be  deposited  in  a layer  thin  compared  with  the  very  narrow 
reaction  zone.  Therefore,  approximating  the  solid-phase  reaction  zone 
as  an  interface  seems  justified.  Figure  1 illustrates  some  of  the 
definitions. 

The  temperature  field  within  the  solid,  occupying  the  region  x<0, 
is  described  by  the  equation 


pc 


3T 

3t 


+ me 


3T  ' 32T  . 3Q 

ax  - 3x2  ax  ' 


CD 


where  T is  the  temperature,  m,  the  mass  burning  rate,  p,  the  solid 
density,  c,  the  solid  heat  capacity  per  unit  mass  and  X,  the  solid 
thermal  conductivity.  The  symbol  Q,  already  defined,  represents  the 
incident  external  radiant  energy  flux  reaching  the  solid  surface,  and 
having  the  spatial  distribution  previously  stated, 

Q(X)  = Qe^  . (2) 

In  general,  Q can  be  a function  of  time  also,  but  this  part  of  the 
study  concerns  only  stimulation  under  steady-state  conditions.  Thus, 
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Figure  1.  Effective  Radiant  Energy  Flux  and  Characteristic 
Lengths  in  Solid. 
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with  Q independent  of  time,  Eq  (1)  becomes 

„ dT  _ , d2T  yx 

dx  x 2 QlJe  • 
dx 


(3) 


The  boundary  conditions  for  Eq  (3)  are  an  interface  condition  involving 
the  flux  A dT/dx  at  x = 0 and 


and 


T = T at  x = 0 , 
s 


T = T.  at  x = -«  > 


(4) 

(5) 


where  Ts  is  the  surface  temperature  and  the  initial  temperature. 

An  interface  energy  conservation  condition,  generally  applicable,  is 

(x  ai)  = m (L  - Ls)  > (6) 

which  in  this  form  has  the  solid-phase  pyrolysis  process  occurring 
within  the  interface.  In  Eq  (6),  L represents  the  sum  of  the  heat  of 
gasification  and  the  heat  conducted  into  the  interior  of  the  solid,  both 
per  unit  mass  of  gasified  solid;  and  Ls  represents  the  heat  of 
gasification  at  temperature  Ts  per  unit  mass  of  gasified  solid. 

In  nondimensional  form,  Eq  (3)  is 


dx 

d£ 


+ e 


o5 


* 


where  x is  the  nondimensional  solid  temperature 

2 2 

(T  - T^mV 
T = QyX  ' 


(7) 


(8a) 


£ the  nondimensional  distance  in  the  solid 
mcx 

6 = ~ ’ 

and  a the  nondimensional  absorption  coefficient  for  the  solid 


(8b) 


_ yX_ 
me 


(8c) 
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The  boundary  conditions  are  now 


t = 0 at  C 1 , 

and  from  Eq  (6)  at  5 = 0, 

dx  <L  - y 
d£  QpX 


The  solution  to  Eq  (7)  with  the  boundary  condition  of  Eq  (9)  is 

r 


t = Ae 


O (1  - a) 


(9) 


(10) 


(ID 


The  boundary  condition  in  Eq  (10)  requires  that 


A = 


m c (L  - Ls) 

QyX 


l 

l-o  * 


(12) 


Therefore,  in  dimensional  form,  the  solution  for  T is 


This  equation  implies,  through  evaluation  at  x = 0 that 


T = T.  + 
s 1 


1 + 2_ 

mc 


Thus,  an  alternative  way  to  write  Eq  (13)  is 


mex 

T = T.  + (T  - T.)  e x + 
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Furthermore,  the  heats  of  gasification  Ls  and  Lso,  with  and  without 
radiant  flux  present,  respectively,  are  related  by  the  thermodynamic 
identity 


L 
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L + (c  - c) 
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(T  - T ) 
v s so' 
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where  Cp  is  the  specific  heat  at  constant  pressure  for  the  gas  and  Ts 
and  Tso  the  surface  temperature  with  and  without  radiant  flux  present, 
respectively.  Figure  2 illustrates  typical  examples  of  radiant  energy 
and  temperature  distributions  in  the  solid  phase. 

C.  Limiting  Cases  of  Behavior 

Eqs  (13)  and  (15)  contain  two  characteristic  lengths,  the  heat 
conduction  length,  given  by  A/mc,  and  the  radiation  absorption  length, 
given  by  1/y.  As  indicated  before,  both  these  lengths  are  assumed 
large  compared  with  the  length  of  the  solid-phase  reaction  zone.  Three 
limiting  cases  occur,  which  cover  the  alternative  situations  that  could 
be  encountered  in  practice. 

1.  Radiation  Absorption  Length  Large  Compared  with  Heat  Conduction 

Length,  (i.e.,  y«  mc/A).  In  this  case  there  are  two  distinct  regions 

predicted  by  Eq  (13)  or  (15).  Restricting  attention  to  the  position  where 

ux 

the  heat  conduction  zone  is  closest  to  the  surface,  the  term  q varies 
here  very  slowly  and  is  nearly  unity.  Therefore,  Eq  (15)  can  be  written 
as 
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This  expression  clearly  indicates  that  the  only  change  produced  by  the 
radiation  is  to  raise  the  effective  initial  temperature  to  the  value 


T 


i eff 


♦ 2-  . 
me 


(18) 


At  the  surface  of  the  solid,  the  temperature  profile  in  the  heat 
conduction  zone,  as  well  as  all  profiles  in  hotter  regions  outside,  are 
the  same  as  they  would  be  if  the  initial  temperature  of  the  propellant 
were  T^  eff. 

2.  Radiation  Absorption  Length  Small  Compared  with  Heat  Conduction 
Length  (i.e.,  y»  mc/A).  In  this  case,  Eq  (13)  simplifies  to  the 
approximate  relation 
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Figure  2.  Calculated  Distributions  of  Radiant  Energy  and 
Temperature  in  Solid  Phase. 
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T = T.  + 
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at  all  positions  except  the  very  narrow  layer  at  the  surface  of  the  solid. 
The  interpretation  here  is  that  instead  of  the  gasification  requiring 
an  amount  of  heat  Ls  at  the  surface,  it  now  requires  a lesser  amount, 

Ls  - (Q/m),  because  of  the  additional  energy  input  at  the  surface.  The 
initial  propellant  temperature  remains  Tj  in  Eq  (19)  and  the  heat 
conduction  zone  maintains  its  single  exponential  structure.  Effectively, 
the  absorption  of  radiation  merely  modifies  the  energetics  at  the 
interface.  Therefore,  provided  the  gas-phase  reaction  is  not  controlling 
and  the  surface  temperature  adjusts  itself,  no  effect  on  burning  rate 
will  occur. 

3.  Radiation  Absorption  Length  Equal  or  Nearly  Equal  to  the  Heat 
Conduction  Length  (i.e.,  y = mc/A  or  y = mc/A) . This  case  separates 
the  two  previous  limiting  cases.  When  the  two  lengths  are  identically 
the  same,  (y  = mc/A),  the  form  of  the  solutions  in  Eqs  (13)  and  (15)  is 
indeterminate.  By  solving  the  problem  through  limiting  procedures,  a 
term  involving  x exp  (mcx/A)  must  appear.  Since  it  is  unlikely  that  the 
two  lengths  will  be  exactly  equal,  the  solution  for  this  case  is  not 
important.  When  the  two  lengths  are  approximately  of  the  same  magnitude, 
the  equations  must  be  used  as  they  stand  without  further  simplification. 

D.  Applicability  of  Equivalence  Principle 

As  indicated  above,  the  condition  y«  mc/A  allows  the  effect  of 
external  radiant  energy  input  to  be  equivalent  in  all  respects  to  an 
increase  in  initial  temperature  by  an  amount  Q/mc,  the  result  given  by 
Eq  (18).  However,  a further  condition  should  be  satisfied  which  concerns 
the  applicability  of  the  steady-state  approximation  itself.  If  the 
propellant  is  too  transparent,  then  this  steady-state  effective  initial 
temperature  might  not  be  achieved  in  the  course  of  an  experiment.  The 
magnitude  of  the  nonsteady  transient  is  derived  from  Eq  (1) . The 
characteristic  time  for  the  solid  to  respond  to  deep  in-depth  radiant 
absorption  can  be  calculated  by  balancing  the  first  and  last  terms 
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in  Eq  (1) . Since  the  characteristic  temperature  rise  that  needs  to  be 
achieved  is  Q/mc,  the  balance  produces  the  characteristic  time 


This  implies  that  if  t*  corresponds  to  the  time  the  external  radiation  is 
applied,  then  achieving  steady-state  behavior  in  the  experiment  requires 
that  y>>  pc/mt*.  Therefore,  the  general  condition  for  the  radiation 
effect  to  be  equivalent  to  an  increase  in  initial  temperature  can  be 
written  as 
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me 
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(21) 


The  upper  condition  on  y given  above  can  be  relaxed  for  many  applications. 
When  y is  comparable  to  mc/A,  the  temperature  profile  in  the  heat 
conduction  zone  of  the  solid  is  altered,  but  the  burning  rate  change  is 
still  equivalent  to  an  initial  temperature  increase,  provided  that  none 
of  the  chemical  reactions  controlling  the  burning  rate  occur  in  this  heat 
conduction  zone.  In  fact,  if  the  chemical  reactions  of  importance  occur 
only  in  the  gas-phase  flame,  even  the  surface  absorption  case,  y»  mc/A, 
is  formally  equivalent  to  an  increase  in  effective  initial  temperature 
insofar  as  the  burning  rate  is  concerned.  The  basic  requirement  for  the 
burning  rate  to  be  described  by  T^  e££  is  that  the  radiant  energy  be 
absorbed  beneath  the  zone  in  which  are  occurring  whatever  chemical 
reactions  that  control  the  burning  rate. 

Whenever  the  radiation  effect  can  be  viewed  as  an  increase  in 
initial  propellant  temperature,  the  situation  possesses  a great  advantage, 
viz.,  there  is  no  requirement  for  a theoretical  combustion  model  capable 
of  predicting  the  burning  rate.  The  effect  of  thermal  radiation  can  be 
obtained  even  if  the  structure  of  the  reaction  zone  is  exceedingly 
complicated.  The  calculation  is  done  through  the  coefficient  of 
temperature  sensitivity,  as  discussed  in  Section  IV.  Of  course,  this 
procedure  implies  that  the  dependence  of  burning  rate  on  initial  temperature 
is  known  experimentally  over  the  pressure  range  of  interest. 
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E.  Burning  Rate  Augmentation 

The  ultimate  effect  of  the  radiant  flux  on  the  combustion  of  the 
propellant  is  to  increase  its  burning  rate,  pressure  conditions 
remaining  constant.  The  magnitude  of  the  burning  rate  augmentation  is 
derived  by  incorporating  the  equivalence  principle  into  the  rate 
expressions.  Since  either  the  gas-phase  or  the  solid-phase  reaction  can 
independently  control  the  burning  rate  under  varying  conditions,  these 
two  alternative  cases  are  examined  separately. 


1.  Control  by  Gas-Phase  Reaction.  The  analysis  of  the  gas-phase 
reaction  zone  is  given  in  Appendix  A.  The  activation  energy  is  assumed 
to  be  high.  Under  gas-phase  control,  even  the  case  of  radiation 
absorption  at  the  solid  surface  satisfies  the  equivalence  rule,  and  can 
be  described  as  an  effect  on  initial  propellant  temperature.  The  analysis 
of  the  gas  phase  is  done  in  such  a way  that  it  can  fit  quasisteadily 
into  a full  time-dependent  analysis  of  the  solid  phase.  The  final 
expression  for  the  gas-phase  flame  model  takes  the  form 

m = CPv/2  e-E/2RTf  * (22) 

where  m is  the  mass  burning  rate,  C,  a constant  independent  of  time,  P, 
the  pressure,  v,  the  order  of  reaction,  E,  the  activation  energy  and  Tf, 
the  flame  temperature.  Eq  (22)  states  that  for  a given  reaction  order, 
the  burning  rate  is  a unique  function  of  pressure  and  flame  temperature. 
The  radiation  will  not  influence  the  combustion  pressure  but  it  will 
change  the  flame  temperature  through  energy  conservation.  Eq  (A8)  of 
Appendix  A can  be  written  dimensionally  as 
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where  H is  the  heat  released  by  the  gas-phase  reaction  per  unit  mass  of 
reactant  consumed;  other  symbols  have  been  previously  defined.  Combining 
Eq  (23)  with  Eqs  (14)  and  (16)  gives  the  overall  energy  conservation 
condition, 
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(24) 
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The  zero  subscript  denotes,  as  before,  the  condition  with  no  external 
radiation  present,  i.e.,  a reference  condition.  Eq  (24)  is  seen  to  be 
identical  with 


T-  * T.  ♦ -9-  , 

f fo  me 

P 

as  expected.  Substitution  of  Eq  (25)  into  Eq  (22)  yields 

E / Q 


(25) 
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Eq  (26)  is  the  theoretical  result  for  burning  rate  augmentation, 
conveniently  expressed  in  terms  of  the  burning  rate  ratio  as  a function 
of  the  radiant  flux  Q and  other  well  defined  parameters.  Most  instances 
of  practical  interest  are  likely  to  have  Q«mc  T^,  hence  Eq  (26)  may 
be  well  approximated  by 


— = exD  l KQ  \ 
mo  \ (m/raoJ  j ’ 


where  the  coefficient  K is  defined  by 
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The  effect  predicted  by  Eq  (27)  is  shown  graphically  in  Figure  3.  This 
plot  can  be  used  to  verify  experimental  results  whenever  the  value  of 
K defined  by  Eq  (28)  is  known. 


For  small  values  of  the  nondimens ional  parameter  KQ,  Eq  (27)  reduces 


to 
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and  for  very  large  values  of  Q,  Eq  (26)  gives  the  limit 
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Figure  3.  Theoretical  Prediction  of  Burning  Rate  Augmentation. 
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Eqs  (27)  and  (29)  indicate  that  for  small  values  of  absorbed  radiant 
flux,  the  augmented  burning  rate,  m,  increases  roughly  linearly  with 
the  stimulus,  in  contrast  with  the  less  than  linear  effect  predicted  at 
progressively  higher  levels  of  stimulation.  At  these  higher  values,  the 
increase  in  burning  rate  causes  the  dependence  on  radiant  flux  to  be  less 
pronounced.  The  derivative  of  Eq  (27)  is 


dm 

dQ 


m K 
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(31) 


which  for  sufficiently  small  values  of  the  parameter  KQ  reduces  to 


dm 

dQ 


= m K . 
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(32) 
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This  result  is  in  agreement  with  the  analysis  of  Coates  and  Kwak,  who 
obtained  reasonable  good  agreement  with  experiment  for  values  of  K that 
seem  realistic. 


2.  Control  by  Solid-Phase  Reaction.  There  are  circumstances  when 
the  reaction  controlling  the  burning  rate  occurs  in  the  solid  phase 
rather  than  in  the  gas  phase.  For  a number  of  double-base  propellants 
this  is  believed  to  occur,  at  least  in  certain  pressure  ranges.  In 
this  case,  the  previous  treatment  is  equally  applicable  with  suitable 
modifications.  A burning  rate  expression  of  the  form 


E./2RT 
= a6  5 


(33) 


takes  the  place  of  Eq  (22) , A being  a constant  that  involves  the 
preexponential  factor,  and  Es  is  the  activation  energy.  The  energy 
feedback  from  any  gas-phase  flame  is  now  assumed  negligible.  Two  other 
conditions  are  that  the  radiation  must  be  absorbed  in  the  propellant 
beyond  the  surface  reaction  zone,  and  that  a certain  amount  of  heat 
must  be  released  by  a solid-phase  reaction.  This  problem  is  extensively 
investigated  in  the  next  Section.  Then  an  adiabatic  energy  conservation 
condition  determines  the  surface  temperature  Ts,  and  the  end  result  is 
an  expression  like  Eq  (33).  Therefore,  Eq  (27)  still  holds,  but  with 
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the  coefficient  K now  defined  as 
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The  new  coefficient  Ks  is  larger  than  the  corresponding  coefficient  X 
for  the  gas  phase,  mainly  because  Tso  is  smaller  than  T£Q,  so  that  a 
stronger  burning  rate  augmentation  can  be  expected  from  the  same 
external  radiant  stimulus  Q under  solid-phase  control. 

Konev'’  worked  out  a different  theory  by  assuming  that  Ts  remains 
Ts0  in  the  presence  of  external  radiation.  This  implies  the  existence 
of  a fixed  gasification  temperature  for  the  solid  surface,  which  seems 
doubtful  because  the  temperature  adjusts  itself  to  maintain  the  steady 
state.  Therefore  his  result  is  questionable.  Of  course,  his  results 
for  the  opposite  case,  deep  in-depth  absorption  of  radiation,  are  not. 
Konev's  later  experiments  appear  to  be  the  most  accurate  of  all  and 
his  reported  data  verify  well  the  predictions  for  the  above  case.  As 
seen  here,  this  is  the  case  of  greatest  interest  for  practical  propellants. 


III.  SOLID-PHASE  REACTION  ZONE 
A.  Inclusion  of  Heat  Generation 

Exothermicity  within  the  solid  phase  may  contribute  significantly 
to  the  burning  rate  particularly  for  double-base  propellants.  Work  by 

collaborators  of  Konev6  has  employed  this  view  to  obtain  a simplified 
description  of  all  aspects  of  the  combustion  of  double-base  propellants, 

including  effects  of  external  radiation.  This  work,  as  well  as  much 
related  work  on  solid-phase  exothermicity,  adopts  a burning  rate 
formula  of  Zeldovich7  which  is  not  exactly  correct  in  the  limit  of  high 
activation  energy.  Because  of  its  interest  and  potential  importance, 
an  improved  treatment  of  solid-phase  reactions  with  absorption  of 
radiation  is  developed  here. 
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For  a solid  in  the  region  x<0,  the  conservation  equations  may  be 
written  in  the  form 


dY 

nw—  = -w  , 
dx  * 


(35) 


where  w is  the  reaction  rate  and  Y,  the  mass  fraction  of  reactant  in 
the  solid,  and 


mc^l  = x ^ + Hw  + Qp  eUX  , 
dx 


(36) 


where  H is  the  heat  released  in  the  solid-phase  reaction  per  unit  mass 
of  reactant  consumed,  other  than  by  activity  of  the  propellant 
stabilizers.  The  remaining  symbols  have  the  same  meaning  as  in  Section 
II.  Reactant  diffusion  within  the  solid  is  neglected,  and  for 
simplicity,  a zero-order  reaction  is  assumed,  allowing  the  reaction  rate 
w to  take  the  form 


-E  /RT 

w = pA  6 


(37) 


where  A is  the  preexponential  factor  in  sec  *,  and  Eg  is  the  activation 
energy  for  the  solid-phase  reaction.  The  boundary  conditions  for 
Eqs  (35)  and  (36)  are 
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The  surface  temperature  is  specified  in  the  traditional  manner,  through 
surface  heat  transfer  and  other  processes  to  be  considered  later. 

Similarly  to  the  gas-phase  analysis,  the  result  seeks  the 
asymptotic  solution  for  large  activation  energy.  There  develops  a 
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convective-diffusive  zone  where  the  solution  to  Eq  (35)  is  Y 
Eq  (36)  reduces  to 

dT  _ d^T 
mc-r-  "“A  o j 
dx  . 2 J 
dx 


1,  and 


(40) 


under  the  assumption  that  no  radiation  is  absorbed  in  the  convective- 
diffusive  zone.  Only  limiting  cases  of  the  analysis  given  in  Section  II 


will  be  used.  The  solution  to  Eq  (40)  is 

mcx 
X 
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T = Ti  + constant  e 


(41) 


and  substituting  for  the  simplest  case,  when  Q = 0 and  the  rate  w is 

localized  at  the  surface,  shows  that  the  constant  is  (T  - T.). 

si 

There  is  reason  to  generalize  Eq  (41)  by  assuming  that  an  exothermic 
reaction,  releasing  a small  amount  of  heat,  occurs  within  this  zone. 

The  reaction  is  interpreted  to  be  the  well  known  reaction  between  the 
stabilizer  and  the  nitroglycerin  or  nitrocellulose. 1,6  The  resulting 
temperature  profile  differs  from  Eq  (41)  and  depends  on  the  rate  of  the 
new  reaction,  whose  activation  energy  is  not  likely  to  be  as  large  as 
Es.  Assuming  that  the  new  reaction  goes  to  completion  early,  only  the 
temperature  gradient  at  the  zone-end  nearest  to  the  propellant  surface  is 
needed.  In  this  region,  if  the  new  reaction  liberates  the  heat  H^,  then 
Eq  (41)  changes  via  energy  conservation  to 
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Eq  (42)  takes  the  form 


T = T„  + (T  - T„)  e 
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The  next  point  to  consider  is  the  way  in  which  radiant  energy  modifies 


(44) 
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Figure  4 
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Solid-Phase  Reaction  Zones  Without  External 
Radiation. 
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Figure  5.  Solid-Phase  Reaction  Zones  With  External 
Radiation. 
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Eq  (44).  Figures  4 and  5 qualitatively  illustrate  the  temperature 
profiles  and  other  definitions  of  the  solid-phase  reaction  zone. 

B.  Limiting  Cases  of  Behavior 

There  are  again  two  final  cases  of  interest.  To  obtain  a general 
expression  capable  of  exhibiting  two  limiting  types  of  behavior,  the 
following  approach  is  utilized.  It  considers  that  out  of  the  total 
radiant  flux  Q reaching  the  surface,  a portion  is  absorbed  in  the 
interior  of  the  solid  as  prescribed  by  the  condition  y«mc/A  and  the 
remainder  Q2,  is  absorbed  as  prescribed  by  the  condition  y>>mc/A,  under 
different  localized  possibilities.  It  will  be  later  ascertained  tha.  he 
radiation  Q2  needs  to  be  treated  in  detail  only  when  its  absorption  takes 
place  in  a layer  of  the  same  or  smaller  length  than  the  surface  reaction 
zone. 


1.  Effect  of  Radiation  Q^,  Having  an  Absorption  Length  Large 
Compared  with  the  Heat  Conduction  Length.  For  y«mc/A,  the  radiation  is 
absorbed  in  a convective-absorptive  zone  deep  inside  the  solid.  Its 
effect  then  appears  in  much  the  same  way  as  Hj.  Therefore,  Eq  (44)  is 
changed  to 
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where  the  constant  of  integration  has  been  evaluated  by  using  the 
requirement  of  Eq  (38b),  that  T = Ts  at  x = 0,  the  lowest-order  matching 
result.  This  case  is  depicted  in  Figure  5 (a). 

2.  Effect  of  Radiation  Q2,  Having  an  Absorption  Length  Small 
Compared  with  the  Heat  Conduction  Length.  For  y>>mc/A,  the  radiation  is 
absorbed  in  a diffusive-absorptive  zone  located  adjacent  to  the 
convective-diffusive  zone,  which  is  illustrated  in  Figure  5 (b) . There 
are  two  limiting  possibilities  depending  upon  whether  the  diffusive- 
absorptive  zone  length  is  large  of  small  compared  with  the  length  of  the 
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reaction  zone,  which  is  located  at  the  solid  surface  in  the  limit  of  high 
activation  energy.  If  the  diffusive-absorptive  zone  length  is  large,  then 
the  reaction  zone  solution  should  be  matched  to  this  instead  of  to  the 
convective-diffusive  zone.  This  case  merits  no  further  analysis,  since 
in  a sense  it  is  intermediate  and  is  not  likely  to  be  encountered  in 
practice.  It  requires  three  distinct  length  scales  to  coexist 
simultaneously  such  that  A/mc  >>  1/p  » A/Bmc,  where  B is  a nondimensional 
activation  energy  defined  in  Appendix  B.  The  remaining  significant 
possibility  is  when  Q2  is  absorbed  in  a length  of  the  same  or  smaller 
size  than  the  reaction  zone  length.  The  full  analysis  is  carried  out  in 
Appendix  B. 


C.  Burning  Rate  Relationship 


The  resulting  burning  rate  expression, 
for  the  complete  solid-phase  reaction  zone, 
sional  form 
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derived  in  Appendix  B, 
takes  the  following  dimen- 
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It  may  be  at  first  surprising  that  Q2  does  not  appear  in  Eq  (46). 

However,  upon  reflection  this  result  is  seen  to  be  entirely  consistent 
with  the  previous  finding  that  for  short  absorption  lengths,  the  radiant 
flux  is  equivalent  to  a decreased  heat  of  gasification.  A total  heat 
balance  would  exhibit  a dependence  of  m upon  Q2,  but  with  the  surface 
temperature  assumed  independent  of  Q2  the  radiation  absorbed  at  the 
surface  zone  is  ineffective  on  burning  rate.  This  result  differs  from 
the  theoretical  result  of  Konev,5  which  contains  an  explicit  relation 
between  m and  Q2.  The  Konev  result  comes  from  an  approximate  treatment 
of  the  reaction  zone,  not  rigorously  valid  in  the  limit  of  high 
activation  energy.  In  particular,  both  the  surface  temperature  and  the 
surface  temperature  gradient  are  specified  in  the  Konev  treatment.  The 
present  analysis  indicates,  (see  Appendix  B),  that  the  surface  temperature 
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gradient  is  determined  by  the  steady-state  solution,  and  that  at  x = 0 


me  \dx  / 
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Essentially,  Konev^  has  not  considered  Eq  (35),  needed  to  ascertain  the 
burning  rate  through  analysis  of  reactant  consumption. 


D.  Generalization  of  the  Zeldovich  Formula 


It  is  of  interest  to  compare  Eq  (46)  with  the  expression  that 
Zeldovich  developed  in  considering  the  standard  case  of  Hj  = 0 and 


Qx  = 0. 


Zeldovich  equation  is' 
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It  is  seen  that  the  RHS  of  Eq  (46)  must  be  multiplied  by  the  term 


h[2c(Ts  - T.)  - h] 


(49) 
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to  recover  the  Zeldovich  equation  for  the  standard  case.  To  recover  the 
generalization  of  the  Zeldovich  formula,  the  RHS  of  Eq  (46)  must  be 
multiplied  by  the  term 
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These  differences  show  that  the  Zeldovich  equation  is  not  rigorously 
valid  in  the  limit  of  high  activation  energy.  In  particular,  it  yields 
in  the  standard  case  finite  burning  rates  for  all  values  of  the  ratio 
c(T  - T.)/H.  Physically,  it  is  clear  that  in  this  case  all  of  these 
treatments  must  become  invalid  when  H ^ c(^s  " T^),  which  corresponds 
to  an  adiabatic  surface  condition  according  to  Eq  (47).  The  general 
outcome  is  that  conditions  for  which 
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are  excluded.  Clearly,  when  the  adiabatic  condition  is  approached,  the 
combustion  wave  will  leave  the  surface  and  propagate  into  the  interior 
of  the  solid,  leaving  behind  essentially  unsteady  conditions  between  the 
wave  and  the  surface.  The  physical  way  to  avoid  this  to  happen  is  to 
allow  the  surface  temperature  Ts  to  increase. 

E.  Overview  of  Theory 

Particularly  for  double-base  propellants,  the  dividing  plane 
between  the  solid  and  gas  phases  is  difficult  to  define.  The  process 
tends  to  be  continuous,  with  the  gas  phase  emerging  clearly  only  towards 
the  end  of  the  fizz  zone.  The  total  amount  of  heat  released  by  the 
end  of  the  fizz  zone  is  estimated  to  be  about  500  cal/g However,  the 
fraction  H of  this  amount  released  in  the  solid  phase  is  difficult  to 
determine.*  For  this  reason,  it  may  be  worth  indicating  that  Eqs  (46) 
and  (47)  can  in  fact  be  applied  for  any  arbitrarily  selected  location  of 
the  solid-gas  interface.  It  is  only  necessary  that  (H  + Hj)  be  the  total 
amount  of  heat  released  per  unit  mass  below  the  selected  interface,  and 
Ql/m  be  the  total  amount  of  radiant  energy  absorbed  per  unit  mass  below 
this  interface.  At  the  end  of  the  fizz  zone,  dT/dx  becomes  very  nearly 
zero,  except  possibly  at  very  high  pressures.  Therefore,  for  certain 
purposes,  it  may  be  desirable  to  select  the  "interface"  x = 0 as  the 
downstream- end  of  the  fizz  zone.  Eq  (47)  then  defines  the  corresponding 
"adiabatic"  flame  temperature  Ts,  and  Eq  (46)  reduces  to  the  classical 

O 

Zeldovich— Frank-Kamenetski  relationship  for  the  flame  speed  . The  gas 
phase  result  of  Appendix  A can  be  viewed  as  a variation  thereof. 

With  this  last  interpretation,  it  is  easy  to  infer  that  the  pressure 
influences  the  reaction  rate;  there  is  no  impediment  to  assigning  A a 


* Reference  5 quotes  a value  of  350  cal/g  for  H and  Reference  6 a value 
of  26  cal/g  for  Hj. 
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pressure  dependence  (typically  A « P1,b  to  obtain  m <r  p0,8)  in  order  to 
obtain  the  pressure  dependence  of  the  burning  rate.  In  the  stricter 
interpretation,  it  is  more  difficult  to  allow  A to  depend  on  pressure, 
because  external  combustion  pressure  should  influence  reaction  rates 
negligibly  in  the  interior  of  the  solid  phase.  Konev^  offers  a plausible 
argument,  keeping  A constant  but  allowing  the  surface  temperature  Ts  to 
vary  with  pressure.  It  is  done  by  maintaining  the  surface  equilibrium 
for  nitroglycerin  in  order  to  recover  the  pressure  dependence  of  burning 
rate.  The  numerical  values  work  out  quite  well,  although  the  required 
heat  of  vaporization  for  nitroglycerin  seems  slightly  high. 

In  so  far  as  the  effect  of  radiant  flux  at  constant  pressure  is 
concerned,  if  the  assumption  of  constant  surface  temperature  is  made^»h, 
as  justified  in  principle  from  the  surface-equilibrium  approach,  then 
it  becomes  clear  from  Eq  (46)  that  the  radiation  Qj  is  equivalent  to  an 
appropriate  increase  in  initial  propellant  temperature,  while  the  radiation 
Q 2 is  totally  ineffective.  The  agreement  with  experiment  in  Reference  6 
supports  this  interpretation  with  Q2  = 0, 

Based  on  the  present  total  study,  this  result  appears  to  have  general 
validity  independent,  for  example,  of  the  assumption  of  constant 
surface  temperature.  It  is  only  necessary  that  some  kinetic  process, 
somewhere  in  either  the  solid  or  gas  phase,  be  responsible  for  determining 
the  burning  rate.  Then  all  radiation  absorbed  below  the  position  where 
this  kinetic  process  occurs  is  equivalent  to  a suitable  increase  in 
initial  propellant  temperature,  while  all  radiation  absorbed  above  this 
position  is  totally  ineffective.  The  variation  from  being  effective  to 
being  ineffective  will  take  place  continuously  and  monotonically. 

IV.  TEMPERATURE  SENSITIVITY 
A.  Significance  of  the  Parameter 

The  variation  of  burning  rate  with  the  initial  temperature  of  the 
propellant  at  constant  pressure  is  represented  by  the  coefficient  of 
temperature  sensitivity.  It  gives  the  dependence  on  temperature,  in 
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contrast  with  the  pressure  exponent  that  gives  the  dependence  on  pressure. 
However,  the  temperature  sensitivity  coefficient  is  amenable  to  formal 
treatment  in  the  context  of  the  theory  of  combustion  thermal  effects. 

Early  studies  on  premixed  gas  flame  theory  produced  the  traditional 
Arrhenius  form  of  the  temperature  sensitivity,  and  the  more  empirical 
forms  which  have  been  found  useful  in  solid  propellant  combustion.1 


When  viewed  as  an  experimental  parameter,  its  customary  role,  the 
temperature  coefficient  offers  the  advantage  of  permitting  calculations 
without  the  need  of  precise  combustion  models.  Of  course,  the  opposite 
applies,  and  the  temperature  sensitivity  cannot  be  predicted  with  any 
degree  of  confidence  unless  accurate  combustion  models  are  known.  If 
this  is  the  case,  the  derivation  is  easily  performed,  as  later  shown  with 
Eq  (46).  By  definition,  the  coefficient  of  temperature  sensitivity, 
denoted  by  (nT)p,  is 


(52a) 

(52b) 


where  subscript  p means  that  the  pressure  is  kept  constant;  r is  the 
regression  rate,  i.e.,  the  usual  measurement  of  solid  propellant  burning 
rate.  From  Eq  (52b),  it  is  obtained  explicitly  that 
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and  introducing  Eq  (18), 
In  r = In 
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Expressed  differently,  Eq  (54)  means  that  the  effect  of  radiant  flux 
can  be  obtained  by  measuring  the  burning  rate  at  elevated  initial 
temperatures.  For  high  radiant  fluxes,  it  may  be  necessary  to  calculate 
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iteratively,  because  the  effective  temperature  depends  on  the  burning 
rate  r in  the  presence  of  external  radiation.  The  new  effective  initial 
temperature  is 


T.  = T.  + -2- 
i xo  rpc 


(55) 


where  is  a reference  initial  temperature,  with  no  external  radiation 
present.  Sometimes  the  following  approximation  is  made: 


r = r 


(Vp 


(Ti  - 


T.  ) 
io' 


(56) 


where  r is  the  burning  rate  at  the  higher  initial  temperature  T^,  and  r , 
the  burning  rate  at  the  initial  temperature  (n^)p  is  the  experimental 

coefficient  for  the  temperature  range  to  that  satisfies  Eq  (56). 
This  equation  is  consistent  with  Eqs  (27)  and  (28)  through  the  relation 
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which  is  the  Arrhenius  form  of  the  temperature  sensitivity  coefficient. 
Changing  Eq  (56)  into 


r = r exp 
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(Vp  7^ 
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gives  a transcendental  expression  in  r,  although  its  inverse. 


rpcln(r/r  ) 

CnT)p  ° > <59) 

is  not  transcendental  in  Q.  An  exact  version  of  Eq  (59),  not  dependent 
on  the  assumption  of  constant  (II^)  is 


Q = rcc  jf  TO"  ' t60) 

o ^ 

When  this  simple  approach  cannot  be  taken,  the  only  recourse  is  to 
delve  into  an  analysis  of  the  reaction  zone. 
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B.  Temperature  and  Pressure  Dependence  on  (fix) p 


In  a general  sense,  the  temperature  sensitivity  of  solid  propellant 
combustion  can  exhibit  dependences  on  both  initial  temperature  and 
pressure,  i.e.,  IIj  * IIj  (T^,  P)  , as  well  as  other  variables.  Double-base 
propellants  display  at  least  two  distinct  types  of  temperature  sensitivity 
behavior,  which  lack  definite  theoretical  explanation.  The  observation 
is  that  certain  analogies  exist  between  pressure  and  temperature 
dependences. 

1.  Linear  Double-Base  Propellants.  The  denomination  "linear"  is 
given  here  to  that  class  of  double-base  propellants  that,  over  an  extended 
pressure  range,  exhibit  a burning  rate  exponent  n and  temperature 
coefficient  (n-jOp  which  remain  essentially  independent  of  mean 
combustion  pressure.  This  is  the  less  complicated  case  and  comprises 
many  well  known  formulations.  Propellants  like  JPN  and  N-4  belong  to 
the  linear  class  under  typical  rocket  operating  pressures  and  initial 
temperatures.  It  is  often  conventional  to  assign  these  propellants  a 
constant  "explosion"  temperature  Te  such  that 


r 


(61) 


where  f(P)  is  a simple  function  of  pressure,  with  constant  coefficient 
b and  exponent  n.  Despite  the  uncomplicated  form  of  Eq  (61),  the 
agreement  with  experimental  data  for  propellants  burning  in  the  linear 
regime  is  remarkably  good.  In  this  case. 
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and  no  pressure  dependence  is  involved.  Figure  6 gives  a plot  of  the 
calculated  values  for  JPN  and  N-4  propellants  as  a function  of  initial 
temperature,  in  terms  of  Eq  (62).  It  is  seen  that  the  temperature 
coefficient  (nT)p  increases  less  rapidly  with  initial  temperature  for  the 
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Figure  6.  Calculated  Temperature  Coefficient  for  JPN  and  N-4 
Double-Base  Propellants. 
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larger  Te.  In  view  of  the  above,  Eq  (61)  can  be  rewritten  as 


r = (nT)p  f(P)  , (63) 

indicating  that  the  burning  rate,  for  the  class  of  propellants  under 
discussion,  increases  linearly  with  the  temperature  coefficient. 


2.  Mesa  and  Plateau  Double-Base  Propellants.  These  formulations 
are  characterized  by  combustion  regimes  where  both  the  pressure  exponent 
and  the  temperature  coefficient  become  strongly  dependent  on  mean  pressure, 
the  latter  retaining  also  a dependence  on  initial  temperature  of  the  type 
discussed  before.  This  is  the  typical  behavior  of  propellants  like  N-5 
and  X-14  which  contain  ingredients,  lead  and  copper  compounds,  responsible 
for  the  mesa  and  plateau  regimes.  This  burning  behavior  is  of  great 
interest  but  no  present  combustion  models  can  fully  interpret  the  complex 
pressure  and  temperature  dependencies  observed  experimentally.  Promising 
chemical  schemes  have  been  recently  proposed  by  Kubota  et  al.^,  which 
give  an  insight  into  the  difficulty  of  the  analytical  part  of  the  problem. 
An  analysis  of  the  affected  reaction  zone  might  lead  to  expressions  of 
the  form 


r = 


-f.cp)  . 

Te(P)-T. 


(64) 


but  solutions  are,  of  course,  unknown.  Restricting  the  phenomenon  to 
the  experimental  observation,  the  following  behavior  appears  to  be 
typical.  For  progressively  higher  pressures,  the  exponent  n first  increases 
and  then  diminishes  rapidly  as  the  mesa  or  plateau  region  is  approached, 
later  recovering  as  the  region  is  surmounted  at  elevated  pressures. 

Along  with  this  variation  of  the  pressure  exponent,  the  temperature 
coefficient  undergoes  an  almost  parallel  drastic  reduction  at  the  mesa 
or  plateau  region;  in  turn,  the  effect  lessens  quickly  thereafter,  as  the 
burning  rate  passes,  this  region  and  increases  at  higher  pressures.  This 
is  the  kind  of  trends  reported  by  Kubota  et  al.9  for  the  N-5  propellant; 
it  is  also  evident  for  the  X-14  propellant.  Using  experimental  burning 
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rate  data  in  conjunction  with  Eq  (56),  the  calculated  (Hj,)p  values 
are  shown  in  Figure  7 as  a function  of  pressure.  Apart  from  the 
theoretical  interest  of  a possible  interrelation  between  pressure 
exponent  and  temperature  coefficient,  the  practical  implications  are 
important.  Propellants  having  small  temperature  coefficients,  whose 
dependence  on  initial  temperature  is  also  low,  are  generally  believed  to 
be  the  most  suitable  for  application.  However,  the  minimization  of  the 
temperature  coefficient  is  bound  to  involve  opposing  requirements,  and 
can  only  be  properly  resolved  for  each  particular  case. 

C.  Effect  of  Stabilizers 

These  compounds  are  widely  used  in  double-base  propellant  formulations 
for  the  purpose  of  preventing  undesirable  chemical  decomposition  during 
storage.  Konev6  and  others  have  advanced  the  idea  that  these  substances 
could  significantly  contribute  to  the  heat  generated  inside  the  solid 
phase  as  the  propellant  bums. 

The  detailed  analysis  of  Section  III  and  Appendix  B was  undertaken 
because  of  the  potential  importance  of  this  effect,  and  its  strong 
relevance  toward  mechanistic  control  of  temperature  sensitivity.  The 
very  comprehensive  nature  of  the  solution  thus  obtained,  principally 
Eq  (46) , permits  a confident  prediction  of  the  influence  of  stabilizers 
on  the  temperature  sensitivity.  Such  a result  is  not  possible  without  a 
solution  like  Eq  (46).  The  overall  indication  is  that  larger  temperature 
sensitivity  (nj)p  is  to  be  expected  as  heat  generation  from  the  stabilizers 
increases.  The  theoretical  result  derived  from  Eq  (46)  is 
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where  y is  a coefficient  given  by 
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This  coefficient  is  determined  by  the  effect  of  the  initial  temperature 
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Figure  7.  Observed  Temperature  Coefficient  for  X-14 
Double-Base  Propellant. 
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on  the  surface  temperature  Tg.  Numerically  y will  vary  between  0 and 
1,  reaching  zero  when  Ts  remains  constant,  although  for  most  propellants 
is  probably  an  intermediate  value.  As  before,  H is  the  heat  generated  in 
the  solid  phase  by  the  main  reaction,  and  Hj  the  heat  generated  by  the 
reaction  pertaining  to  the  stabilizers.  Since  the  surface  temperature  of 
burning  double-base  propellants  is  not  that  high,  about  550°K,  Eq  (65) 
shows  that  the  temperature  coefficient  can  rapidly  rise  from  its  Arrhenius 
datum  level  even  for  modest  increases  of  either  of  the  two  solid-phase 
heat  generation  terms,  but  primarily  by  increases  from  the  stabilizers' 
term.  The  effect  of  increments  is  twice  as  strong.  Eq  (65)  equally 
shows  the  influence  of  the  other  variables  of  interest  and  suggests 
plausible  ways  to  reduce  the  temperature  sensitivity  coefficient  when 
the  burning  rate  is  under  solid-phase  reaction  control. 

In  view  of  the  practical  significance  of  this  result,  experiments 
ought  to  be  done  to  confirm  the  theoretical  prediction.  The  experiment 
should  be  as  follows:  (1)  Selection  of  a linear  propellant  (e.g„,  JPN) 

formulated  with  zero,  low,  medium  and  high  concentrations  of  stabilizer; 
(2)  determination  of  the  heat  generation  terms  from  thermograms  run  at 
fast  temperature- time  conditions;  (3)  determination  of  r(T^,P)  for  each 
concentration  of  stabilizer,  in  order  to  calculate  the  temperature 
coefficients;  and  (4)  comparison  of  results  on  the  basis  of  stabilizer 
concentration  for  each  combination  of  pressure  and  initial  temperature. 

The  measured  heat  generation  terms  might  well  not  be  the  same  as  during 
combustion  but  should  give  an  estimation  of  relative  magnitude.  Also, 
measurement  of  absorption  coefficients  would  resolve  the  effect  that 
stabilizer  concentration  might  possibly  have  on  this  parameter. 
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V.  COMPARISON  WITH  EXPERIMENT 


A.  Calculation  Procedure 

In  view  of  the  theoretical  results,  the  best  way  to  confront 
experimental  evidence  is  on  the  basis  that  the  radiant  flux  is  equivalent 
to  an  increase  in  initial  propellant  temperature.  The  experimental 
verification  of  Eq  (27)  is  not  straightforward  because,  for  a given 
propellant,  there  is  no  a priori  knowledge  on  the  exact  value  of  the 
parameters  K or  Kg,  especially  for  a wide  pressure  range.  Therefore, 
the  general  form  of  the  relationship  that  should  be  tested  is  as  follows. 

The  experimentally-controllable  variables  influencing  the  burning 
rate  r are  the  pressure,  p,  the  initial  temperature,  T^,  and  the 
radiant  flux  absorbed,  Q.  This  dependence  can  be  indicated  by  the 
functional  relationship  r = r (P,  T^,  Q) . The  equivalence  of  Q to  an 
increase  in  initial  temperature  is  expressed  by  the  functional  relation- 
ship 

r (P,  T.,  Q)  * r (P,  T.  + ^ , 0).  (67) 

Calculations  are  performed  most  conveniently  by  first  selecting  a value 

of  pressure.  For  this  pressure,  and  any  other  selected  pressure 

thereafter,  a tabulation  of  r as  a function  of  L at  Q = 0,  and  a 

tabulation  of  r as  a function  of  Q at  a fixed  reference  T^,  say 

Tj.  = T^q,  must  be  available.  In  this  manner,  the  effect  of  Q can  be 

predicted  from  the  r(P,  T^)  data  and  then  compared  with  the  effect 

observed  experimentally  at  the  reference  temperature  T^q.  Thus,  the 

objective  is  to  find  r(P,  T^q,  Q) . From  the  tabulation  of  r as  a 

function  of  T. , including  T.  , r is  known  as  a function  of  T.  - T.  , 
i’  B io'  l io’ 

and  then  a value  of  Q is  calculated  for  each  entry  from  the  equation 

Q = rpc  (T.  - Tio)  . (68) 

This  enables  the  derived  data  of  r as  a function  of  Q,  at  each  pressure, 
to  be  generated. 
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As  explained  before,  the  most  suitable  double-base  propellants 
for  this  procedure  are  those  exhibiting  a linear  burning  regime  at  the 
pressures  of  interest,  and  whose  temperature  sensitivity  coefficients 
are  assigned  a specific  form,  for  example,  the  one  given  by  Eq  (62). 
Therefore,  the  relation  of  burning  rate  as  a function  of  initial 
temperature  can  be  used  to  calculate  the  effect  of  radiant  flux  on 
burning  rate  for  propellants  like  JPN  and  N-4,  for  a wide  range  of 
pressures.  Explicit  predictions  can  be  obtained  for  this  class  of 
propellants  by  noting  that  Eq  (60)  reduces  in  this  case  to 
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for  fixed  pressure  conditions,  and  to 
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for  variable  pressure.  These  equations  predict  that  the  burning  rate 
increases  linearly  with  radiant  flux,  and  that  sizeable  radiant  fluxes 
are  needed  to  increase  the  burning  rate  substantially,  particularly  at 
high  pressures. 


B.  Experimental  Measurements 

Some  test  results  have  been  obtained  with  the  apparatus  depicted 
in  Figure  8.  It  consists  of  a source  of  radiant  energy,  a constant 
pressure  combustion  chamber,  and  optical  elements  to  direct,  divert 
and  refocus  the  radiation  beam  on  to  the  propellant  surface.  The  chamber 
operates  essentially  as  a strand-burner.  Windows  have  been  added  to 
permit  entrance  of  the  radiation  beam  and  for  photography. 

The  radiant  flux  is  measured  at  the  position  normally  occupied  by 
the  propellant  sample  during  initiation  of  burning.  The  measurement  is 
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Figure  8.  Experimental  Apparatus. 


made  by  using  a water-cooled  "Thermogage"  continuous  output  calorimeter. 
No  correction  for  gas  absorption  has  been  attempted  and  the  nominal 
Q value  so  measured  is  carried  into  the  calculations  given  later.  The 
initial  temperature  of  the  propellant  sample  is  monitored  after 
pressurization  of  the  chamber  with  a thermocouple  placed  adjacent  to  its 
surface.  The  propellant  sample  is  a cylinder  1.2  cm  in  diameter  and 
2.4  cm  in  length,  and  otherwise  prepared  and  mounted  inside  the  chamber 
much  in  the  same  way  as  it  is  customary  in  strand  burner  operation. 

The  burning  of  the  propellant  is  recorded  photographically  by  means  of 
a high  speed  "Hycam"  camera,  outside  illumination  being  provided  by  a 
mercury  light.  The  burning  rates  with  and  without  external  radiation 
present  are  determined  from  timed  readings  of  the  regressing  surface, 
taken  by  inspection  of  the  photographic  record.  The  radiant  flux  is  not 
monitored  during  the  course  of  the  run;  the  diameter  of  the  beam  is 
initially  about  3mm  but  some  defocusing  must  occur  as  the  surface 
regresses . 

Many  difficulties  were  encountered  in  performing  the  experiments, 
and  the  number  of  usable  points  eventually  obtained  is  limited  to  a 
rather  low  Q value.  Prevalent  among  the  various  difficulties  is  the 
requirement  to  control  the  initial  temperature  of  the  propellant 
accurately  once  a reference  temperature  has  been  selected.  This  is 
particularly  critical  when  stimulating  at  low  levels  of  radiation.  It 
is  clear  that  even  small  variations  of  initial  temperature  can  have 
effects  comparable  to  that  of  the  external  radiant  flux,  when  only  small 
stimulations  are  available.  Such  a condition  is  a direct  consequence  of 
Eq  (68) , which  is  shown  in  Figure  9.  The  technique  for  measuring  the 
heat  feedback  from  the  gas  flame  alone  was  also  intensively  investigated 
at  the  beginning  of  the  program.  Here,  the  uncertain  sensor  response 
and  light  guide  arrangement  precluded  the  acquisition  of  reliable  data. 

C.  Verification  of  Results 

The  original  experimental  plan  included  the  testing  of  several 
propellants  over  a wide  range  of  combustion  pressures  and  stimulations 
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Figure  9.  Equivalent  Increase  in  Initial  Propellant  Temperature 
for  Various  Burning  Rates. 


47 


of  radiant  flux.  Only  the  JPN  propellant  could  be  tested,  at  a rather 
low  flux  level,  but  the  data  points  extend  to  higher  pressures  than 
any  previous  results.  These  are  listed  in  Table  1,  which  also  contains 
the  information  on  burning  rate  and  propellant  physical  properties  needed 
to  calculate  the  predicted  results.  The  burning  rate  determinations 
given  in  Table  1 are  values  derived  over  many  readings  of  the  film  record. 
Additional  data  points  were  obtained  at  atmospheric  pressure  as  a 
function  of  radiant  flux  stimulation.  The  results  are  presented  in  terms 
of  the  ratio  r/r  , where  rQ  is  the  normal  burning  rate  and  r the 
radiation-enhanced  burning  rate  both  at  the  same  pressure. 

The  experimental  results  for  the  JPN  propellant  are  shown  in 
Figures  10  and  11,  together  with  the  predictions  calculated  from  Eqs  (70) 
and  (71)  respectively.  Due  to  the  intrinsic  difficulty  of  distinguishing 
such  small  variations  in  burning  rate,  of  the  order  of  the  precision 
of  the  experiment  itself,  the  observed  data  scatter  is  not  surprising. 

The  data  reported  in  Reference  6 were  also  analyzed  and  excellent 
agreement  between  observed  and  calculated  results  was  found.  On  the 
basis  of  all  available  data,  the  view  that,  in  so  far  as  its  effect  on 
burning  rate  is  concerned,  radiant  flux  does  nothing  but  effectively 
elevate  the  initial  temperature  of  the  propellant,  appears  to  be 
supported  quite  well. 
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Figure  10.  Observed  and  Calculated  Burning  Rate  Augmentation 
for  JPN  Propellant  at  Atmospheric  Pressure. 
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Figure  11.  Observed  and  Calculated  Burning  Rate  Augmentation 
for  JPN  Propellant  at  Elevated  Pressures. 
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Table  I 


Burning  Rate  Augmentation  for  JPN  Propellant 


Q = 0 

Q « 1.25  x 105 

CM 

a 

3: 

0.79  x 106  N/m2 

0.79  x 106 

N/m2 

Test 

# 

Init  Temp. 
T.,  °K 

Bum  Rate 
r , cm/sec 

Test 

# 

Init.  Temp 

Bum  Rate 
r,  cm/sec 

51 

295 

0.439 

54 

294 

0.409 

52 

296 

0.318 

55 

294 

0.400 

53 

296 

0.368 

56 

296 

0.391 

1.48  x 106  N/m2 

1.48  x 106 

CM 

a 

40 

294 

0.465 

49 

294 

0.587 

41 

293 

0.505 

50 

295 

0.485 

63 

295 

0.470 

2.86  x 106  N/m2 

2.86  x 106 

N/m2 

37 

295 

0.864 

44 

294 

0.932 

42 

292 

0.724 

47 

293 

0.813 

67 

295 

0.818 

48 

295 

0.787 

65 

295 

0.800 

5.62  x 106  N/m2 

5.62  x 106 

N/m2 

58 

296 

1.24 

60 

295 

1.22 

59 

296 

1.30 

61 

295 

1.26 

62 

295 

1.24 

Pressure 
P,  N/m^ 

Average  Results 

cm/ sec 

r 

cm/sec 

r/ro 

0.79  x 106 

0.375 

0.400 

1.067 

1.48  x 106 

0.485 

0.514 

1.059 

2.86  x 106 

0.802 

0.833 

1.038 

5.62  x 106 

1.27 

1.24 

0.98 
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Table  I. 


Burning  Rate  Augmentation  for  JPN  Propellant  (Cont'd) 


Data  for  Prediction  of  Results 


r 


5.97  x 10'5  0.70 

535  - T. 

l 


m/sec 


(P,  N/m2;  T.,  K) 


p = 1.6  x 103  Kg/m3 
c = 1.46  x 103  J/Kg  K 
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VI.  DISCUSSION  AND  CONCLUSIONS 


The  analytical  model  developed  here  predicts  realistically  the 
effect  of  external  thermal  radiation  on  the  burning  rate  of  double- 
base propellants.  Except  for  noted  instances,  the  traditional 
hypothesis,  that  the  effect  can  be  treated  as  an  increase  in  initial 
propellant  temperature,  has  been  found  generally  applicable.  Therefore, 
the  hypothesis  is  more  properly  referred  to  as  an  equivalence  principle. 

The  experimental  evidence  shows  satisfactory  agreement  between 
observed  and  predicted  results,  supporting  the  theoretical  conclusion 
that  the  effect  is  calorific  in  nature  and  not  directly  dependent  upon 
the  chemistry  of  the  combustion  process.  Owing  to  experimental 
difficulties,  the  amount  of  data  collected  in  the  BRL  facility  was  small. 
The  need  for  extensive  new  data,  however,  is  not  critical  because  of 
available  results  from  other  sources.  Several  other  manifestations  of 
thermal  effects  have  been  investigated,  in  particular  the  role  of 
stabilizers  in  temperature  sensitivity.  Experiments  to  elucidate  this 
important  question  are  being  planned. 

The  main  conclusion  is  that  the  radiant  energy  absorbed  by  the 
propellant  has  the  effect  of  raising  the  effective  combustion  temperature 
and  therefore  of  enhancing  the  burning  rate  in  a manner  equivalent  to  an 
increase  in  initial  temperature.  For  typical  double-base  propellants, 
the  radiation  effect  will  remain  relatively  constant  with  time;  in  rocket 
chambers  the  increase  in  burning  rate  over  the  strand-burner  rate  is 
normally  less  than  5%.  On  the  contrary,  transparent  propellants  will 
sustain  a progressive  effect  on  burning  rate;  if  large  amounts  of  radiant 
energy  are  absorbed  deep  in  the  interior  of  the  solid,  other  effects 
like  cracking  and  fissuring  can  appear.  Such  effects  are  not  included 
here. 

The  conclusions  are  also  far  reaching  in  regard  to  the  effect  of 
stabilizers.  The  rise  in  temperature  by  heat  conduction  within  the  solid 
initiates  the  decomposition  of  nitroglycerin  and  nitrocellulose  according 
to  traditional  notions.  These  are  exothermic  processes  liberating  an 
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amount  of  heat  represented  in  the  theory  by  the  term  H.  A second  type 
of  reaction  within  the  solid  that  can  contribute  equally  substantial 
amounts  of  heat  is  the  one  postulated  between  the  nitrate  ester  decomposition 
products  and  the  stabilizers,  represented  in  the  theory  by  the  term  Hj. 

The  function  of  these  compounds  is  well  known  and  has  been  mentioned 
before.  The  fact  is  that  self  heating  must  increase  with  stabilizer 
concentration. 

The  problem  of  temperature  sensitivity  is  basically  part  of  the 
general  theory  of  thermal  effects.  The  concern  here  has  been  to 
emphasize  the  role  of  this  parameter  in  the  theory  and  to  show  that  it 
can  be  predicted  with  confidence  only  to  the  extent  that  the  combustion 
model  used  is  realistic.  The  scope  of  the  analysis  has  lead  to  a 
generalization  of  the  Zeldovich7  formula  and  temperature  sensitivity 
prediction  in  which  the  effect  of  stabilizer  concentration,  among  other 
parameters,  is  clearly  indicated.  Increasing  the  concentration  of  a 
stabilizer  that  reacts  at  low  temperature  increases  the  temperature 
sensitivity  through  the  term  H^;  however,  twice  as  large  of  an  increase 
in  the  heat  release  H of  the  main  solid  reaction  is  needed  to  produce  an 
equivalent  augmentation.  Therefore,  stabilizer  modifications  are 
predicted  to  have  a large  effect. 

Future  analytical  studies  will  consider  the  influence  of  radiant 
flux  when  such  a stimulus  is  applied  in  the  oscillatory  mode. 
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APPENDIX  A 

ANALYSIS  OF  GAS-PHASE  REACTION  ZONE 


Using  the  present  notation,  the  gas-phase  flame  exists  in  the  region 
x>0.  The  quasisteady  conservation  equations  for  energy  and  for  reactant 
mass  fraction  are 


me  4 * « HYP*  e-E/RT 

P dx  g dx2  g 


and 


dY 


d2Y 


m -r—  = p D _ 
dx  g dx2 


- BYPV  © 


•E/RT 


(Al) 


(A2) 


respectively,  where  m is  the  mass  burning  rate,  Cp,  the  heat  capacity  at 

constant  pressure,  X , the  thermal  conductivity,  p , the  density,  H , 

S o E 

the  heat  released  per  unit  mass  of  reactant  consumed,  B,  the  preexponential 
factor  in  the  Arrhenius  rate  expression,  Y,  the  reactant  mass  fraction, 
v,  the  order  of  reaction,  E,  the  activation  energy  and  D,  the  diffusion 
coefficient.  The  boundary  conditions  for  Eqs  (Al)  and  (A2)  are 


T = T at 
s 

X = 

0 

(A3a) 

m(l  - Y)  = 

V 

4^-  at  x = 0 
dx 

(A3b) 

r a dT 

mL  = X -r~ 
g dx 

at 

x = 0 , 

(A3c) 

Y = 0 at  x 

= CO 

J 

(A4) 

where  L in  Eq  (A3c)  is  the  sum  of  the  heat  of  gasification  and  the  heat 
conducted  into  the  interior  of  the  solid,  both  per  unit  mass  of  reactant 
consumed.  A number  of  simplifying  assumptions  are  introduced;  for 
example,  the  thermal  conductivity  and  the  diffusivity  p D have  been 
assumed  constant.  The  reaction  has  been  assumed  to  be  first  order  with 
respect  to  the  reactant,  although  an  arbitrary  order  v has  been  permitted 
insofar  as  the  pressure  dependence  is  concerned.  It  is  later  introduced 
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that  pD=A/c  = a is  a constant,  justified  since  the  Lewis  number  is 
8 8 P ai 

assumed  to  be  unity  . The  dependent  variables,  as  well  as  the  pressure, 
P,  the  mass  rate,  m,  the  surface  temperature,  Tg  and  the  heat,  L,  are  all 
considered  to  be  functions  of  time;  the  heat  transfer  to  the  solid  enters 
through  L. 

Dividing  Eq  (Al)  by  H and  adding  the  result  to  Eq  (A2)  gives 

o 


d(Y  + 0)  d (Y  + 0) 

dx  ,2 

dx 


(A5) 


where  0 = c T/H  is  a nondimensional  gas  temperature.  The  only  solution 
P S 

to  Eq  (A5)  remaining  bounded  at  infinity  is  a constant.  Evaluating 
the  constant  by  using  Eq  (A4)  gives 


Y = 0£  - 0 , 


(A6) 


where  the  time  dependent  flame  temperature  T£  is  nondimensionalized  as 
0£  = CpT£/Hg.  Substitution  of  Eq  (A6)  into  Eq  (A3b)  reveals  that  at 
x = 0, 


m(l  +0  - 0 _)  = a » 

s f dx  * 


(A7) 


where  ©s  applies  to  the  surface.  Eqs  (A3c)  and  (A7)  lead  to  the  solution 


0,,  = 1 + 0 - z 

r s 


(A8) 


where  z = L/H  is  the  nondimensional  parameter  accounting  for  the  heat 

o 

flux  into  the  solid.  Substitution  of  Eq  (A6)  into  Eq  (Al)  produces  the 
equation 


" a! = “ 7§  * BpV  (®f  - e>  e‘EVe 

dx 


(A9) 


where  E'  = Ec  /RH  is  a nondimensional  activation  energy.  The  boundary 
P S 

conditions  for  Eq  (A9)  are 


0 = 0 at  x = 0 , 
s 


(AlOa) 
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from  Eq  (A3c) , 

d0  ^ 

a * mz  at  x = 0 , 

(AlOb) 

and 

0 = 0£  at  x = “ . 

(All) 

A2 

From  the  analysis  of  Bush  and  Fendell  , it  is  known  that  in  the  limit 
of  large  activation  energy,  there  is  a convective-diffusive  zone 
occupying  a position  0 which  is  not  very  close  to  0^.  See  Figure  A1 
for  definitions.  In  the  convective-diffusive  zone,  the  last  term  in  Eq 
(A9)  is  negligible,  and  the  general  solution  is 

0 = a + b e^01  , (A12) 

where  a and  b are  constants  of  integration.  The  boundary  conditions  of 
Eqs  (AlOa)  and  (AlOb)  may  be  used  to  evaluate  the  constants  of 
integration,  giving 

0 = 0s  + z (e^0  - 1)  . (A13) 

In  this  solution,  it  is  seen  by  using  Eq  (A8)  that  0 = 0^  at  x = x^, 
where 

xf = (;)  ln  (i)  • (A14) 

is  to  lowest  order  (E*  -*■  «)  the  distance  between  the  solid  surface  and 
the  gas-phase  reaction  zone.  Differentiating  Eq  (A13)  and  using  Eq 
(A14)  gives 

a!"  I «*-*£•  <A15> 

For  large  activation  energy,  there  exists  in  lowest  order  a 
reactive-diffusive  zone  about  the  position  X£.  In  this  zone,  the  first 
term  in  Eq  (A9)  is  negligible.  The  appropriate  stretched  spatial 
variable  is 

Bm(x  - xf) 

n = T-*-  , (A16) 
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(A17) 


where  the  large  expansion  parameter  is 

EH 

9 


B = ILL  s 
B “ 2 2 
9f  cpRTf 


and  Eq  (A9)  reduces  in  lowest  order  to 


dn2 


Aye"y  , 


where 


y = 6 (0f  - 0) 


is  the  stretched  nondimensional  temperature,  and 

7©f 


• ■ (#K" 


(A18) 


(A19) 


(A20) 


The  boundary  conditions  for  Eq  (18) , derived  by  matching  from  Eq  (A15) 
and  from  the  condition  of  Eq  (All),  are 


y -+  0 as 


(A21) 


and 


d£ 

dn 


as 


(A22) 


Eqs  (A18) , (A21)  and  (A22)  define  the  classical  problem  for  the  flame- 

speed  eigenvalue  for  a premixed  laminar  flame.  The  solution  is 
A?  A3 

A = 1/2  * . Therefore  the  result  from  Eq  (A20)  is 


m = 


2aR2T^c2BPV 

LE 

2 2 
EH 

g 

„2 


1/2 


-E/2RT, 


(A23) 


Although  a factor  T £ emerges  from  the  square  root,  it  can  be  eliminated 
by  assigning  suitable  temperature  dependences  to  a or  B;  in  any  event, 
for  large  activation  energy  the  effect  of  that  factor  is  small  compared 
with  the  effect  of  the  exponential  term.  Hence,  Eq  (A23)  can  be  written 
as 

m = CPv/2e"E/2RTf  , (A24) 
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where  the  coefficient  C is  independent  of  time.  Eq  (A24)  is  the  burning 
rate  prediction  for  the  gas-phase  flame. 
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APPENDIX  B 

ANALYSIS  OF  SOLID-PHASE  REACTION  ZONE 


Asymptotic  considerations  show  that  a reactive-diffusive  zone 
exists  at  the  surface  of  the  solid  phase  when  the  activation  energy  is 
large.  The  definitions  are  the  same  as  in  Figures  4 and  5.  In  this 
reactive-diffusive  zone,  the  proper  spatial  variable  is  defined  by 

(Bl) 


„ _ Bmcx 
n = * — • 


where  in  turn,  the  nondimensional  activation  energy  6 is  defined  by 
E 


B h 


RT 


(B2) 


In  this  zone,  the  relevant  temperature  variable  £ is  defined  by 
B(Ts  - T) 


(B3) 


and  the  following  nondimensional  groups  are  introduced  for  transformation 
of  the  differential  equations: 


and 


H 

cT 

(B4) 

s 

yXRT 

s 

mcEs  * 

(B5) 

yXRQ. 

L 

2 2_ 

(B6) 

Expressing  the  nondimensional  burning  rate  in  the  form 

1/2 


Y = m 


E c 
s 

-E  /RT 

RT  XpAe  S S 
s 


(B7) 


Eq  (35), 


dY  4 -VRT 
m = -pA  e 
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becomes 


dY  _ 1_  p 

dn  ' t2  c 


Q-s/b) 


and  Eq  (36) , 


mc  S = x ~t  + Hw  + Qv  eMX 


becomes 


d\  1 d£  = _h_  e(1"C/P) 

dn2  e ~ *2 


+ be 


-an 


The  boundary  conditions  for  Eqs  (B8)  and  (B9)  are 
Y = 0 at  n = 0 , 


and 


? = 0 at  n = 0 . 


„-l 


(B8) 


(B9) 


(BlOa) 

(BlOb) 


To  lowest  order  in  the  small  parameter  8 , matching  conditions  are 

found  to  be 


and 


Y ->  1 as  n *► 


d? 


-r-2-  C as  n -»•  00  . 
dn 

Here  the  value  of  the  constant  C is  derived  from  Eq  (45)  with  the 
result  that 


(x  -T.-i.5i) 

c = Vs  i c me  / 


(Blla) 

(Bllb) 


(B12) 


Since  solutions  are  sought  only  to  lowest  order  in  the  small  parameter 
B Eq  (8B)  simplifies  to 


dY  _ 1 

dn  ~ ,2  e 

and  Eq  (B9)  simplifies  to 


(B13) 
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(B14) 


Ji_  e-c  + be'an 


The  integral  of  Eq  (B13),  subject  to  the  previously  specified  boundary 
conditions  Eqs  (B10a,b)  and  Eqs  (Blla,b),  provides  an  expression  for  the 
burning  rate. 


2 


dn  . 


(B15) 


The  remaining  problem  is  to  solve  Eq  (B14),  subject  to  the  previously 

specified  boundary  conditions  of  Eqs  (BlOb)  and  (Bllb) . With  the 

_an 

absorptive  term,  b e , included  in  full,  Eq  (B14)  cannot  be  solved 
analytically.  Therefore,  only  the  limiting  case  a.  •*■  <*>  will  be  considered, 
as  explained  in  Section  III.  This  produces  a thin  surface  layer  in  which 
all  the  radiation  is  absorbed,  while  ? changes  very  little.  Within  this 
layer,  the  appropriate  spatial  variable  is 

e * an  , (B16) 

and  the  appropriate  temperature  variable  is 


Thus  Eq  (B14)  becomes 

=*  _h  e -b<f>/a2  + e-e  . (B18) 

de^  'Tb 


As  a + «,  the  last  term  of  Eq  (B18)  tends  towards  a constant,  and  the 
solution  becomes 


♦ = 


2 

E 


(1  - ©G)  + constant  e. 


(B19) 


The  constant  is  to  be  determined  by  matching  with  the  reactive-diffusive 
zone.  It  is  seen  from  Eq  (B18)  that,  to  lowest  order  in  a”  , the 
structure  of  the  reactive-diffusive  zone  is  unaffected  by  the  presence 
of  the  absorption  layer. 
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For  the  reactive-diffusive  zone,  Eq  (B14)  is 


(B20) 


with  lowest  order  boundary  conditions  specified  again  by  Eqs  (BlOb)  and 
(Bllb) , The  first  integral  of  Eq  (B20)  is 


1 

2 


-e 

e 


2 


(B21) 


The  integral  of  Eq  (B21)  is 


n = 


cc2-|eV/2  ' 

¥ 


From  Eq  (B15)  it  follows  that 


(B22) 


whence 


dC 


e~C  dc 

CC2.lhe^l/2 


(B23) 


CB24) 


This  is  the  burning  rate  expression,  which  by  means  of  Eqs  (B7)  and  (B12), 
can  be  written  dimensionally  as  "* 


2 


m 


E 

s 


RT2XpAe_Es/RTs 
s 

feds  - V - £ - Hj 


(B25) 
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This  analysis  shows  that  the  surface  temperature  gradient  is  determined 
by  the  solution,  because  the  derivative  of  Eq  (B22)  implies  that  at  x = 0 


so  that  Eq  (47)  follows.  The  significance  of  all  these  theoretical 
developments  is  presented  in  Sections  III  and  IV. 
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